Line defects in Three dimensional Symmetry Protected Topological Phases 
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A 3d symmetry protected topological phase, by definition must have symmetry protected non- 
trivial boundary states, namely its 2d boundary must be either gapless or degenerate. In this work 
we demonstrate that once we couple a 3d SPT phase to a lattice dynamical Z2 gauge field, in many 
cases the Z2 vison loop excitation (line defect) can be viewed as a "Id boundary" of the 3d SPT 
phase, and this line defect is guaranteed to have gapless or degenerate spectrum, which is also 
protected by the symmetry of the SPT phase. 

PACS numbers: 



In the last few years, motivated by the discovery of 
free fermion topological insulators protected by time- 
reversal symmetry [iHSl, a new class of quantum dis- 
ordered states, the so called symmetry protected topo- 
logical (SPT) phases was proposed T!, 8||. Unlike intrin- 
sic topological phases such as fractional quantum Hall 
states, a SPT phase is only nontrivial when the system 
has certain symmetry G. A d— dimensional SPT phase 
must have a fully gapped and nondegenerate spectrum 
in the bulk, and also a gapless or degenerate spectrum 
on its d — 1 dimensional boundary, when and only when 
the Hamiltonian of the system (both in the bulk and the 
boundary) has symmetry G. In the last two years, SPT 
phase has emerged as a new subfield of condensed mat- 
ter theory, and it has attracted a lot of attentions and 
efforts [JliH. 

Based on the definition of SPT phases, the 2d bound- 
ary of a 3d SPT phase must have nontrivial spectrum. 
But the properties of a Id boundary (or Id hue defect) 
in a 3d SPT has not been studied yet. Line defects in or- 
dinary topological insulators have been discussed before, 
and it was pointed out that these line defects do carry 
gapless modes localized along the defects 22|, |23| . In this 
work we will study one type of line defects in strongly 
interacting 3d bosonic SPT phases, and we will conclude 
that in many cases, this line defect in a 3d SPT phase 
does lead to gapless or degenerate spectrum. 

Since so far we do not have explicit lattice model for 
most of the SPT phases under study, our work will be 
based on the effective field theory of SPT phases. Trivial 
quantum disordered phases can be described as the dis- 
ordered phase of either a Ginzburg-Landau field theory, 
or a semiclassical nonlinear sigma model (NLSM) defined 
with an order parameter. SPT phases have the same bulk 
spectrum and bulk phase diagram as a trivial system, so 
they can still be described by NLSMs, and their nontriv- 
ial boundary spectrum can be captured with a topolog- 
ical 8— term in the bulk [iSJ, [ij]. It was demonstrated 
that the NLSM plus an appropriate topological 0— term 
not only leads to nontrivial boundary physics [24| , it also 
gives us the correct ground state wave function of the 
SPT phase [l^| . In this work wc will focus on several 3d 



SPT phases that are described by the same effective field 
theory, which is a 0(5) Nonlinear Sigma model with a 
0— term at 8 = 2tt: 



5= / d^xdT -{df^nf 
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Here n is a five-component order parameter with unit 
length. Although these SPT phases share the same effec- 
tive field theory, the vector order parameter n transforms 
differently under symmetry in different SPT classes. 

The order parameter n corresponds to certain opera- 
tors on the lattice scales, such as spin and boson oper- 
ators. As long as the symmetry of the SPT phase con- 
tains a Z2 center, i.e. a Z2 subgroup that commutes with 
all the other group elements, we can always modify the 
Hamiltonian by coupling the lattice operators to a dy- 
namical Z2 gauge field on the lattice. Since the matter 
field n is disordered and gapped in the bulk of these SPT 
phases, the Z2 gauge field can have a deconfined phase, 
which is the phase we will focus on in this paper. The de- 
confined phase of a Z2 gauge field introduces line defect 
in the system, which is the vison loop of the Z2 gauge 
field, i.e. a tt— flux loop in the system. We will argue that 
in many SPT phases described by Eq. [Jl the vison loop 
has a nontrivial spectrum, namely it is either gapless or 
degenerate. 

Example 1: 3d SPT with [U{1) x U{1)] x Zj 
Let us start with a simple example of 3d SPT phase 
with [[/(I) X C/(l)] X Z2 symmetry, where Zj is the time- 
reversal symmetry. This SPT phase is described by Eq. [1] 



hi and n 



62 are two independent 



boson fields, and iT' — <j) is_an Ising order parameter that 
changes sign under Zj [13( : 



Zi : 61 ^ 61, 



62 -^ 62 



-> 



(2) 



Eq. [T]has an enlarged SO (5) symmetry, but we can turn 
on extra terms in Eq.[T] which reduce this symmetry down 
to physical symmetry [C/(l) x U{1)] x Z2 . 

We will focus on the SPT phase, namely the phase 
where the five component order parameter n is com- 
pletely disordered. We can couple bi to a Z2 gauge field. 




FIG. 1: A vison loop in example 1 and 2 is bound with a half 
vortex line of 6 ~ n^+in^, which leads to a 1+ld 0(3) NLSM 
with O = TT along the vison loop. 



and let us assume this Z2 gauge field is deep in its de- 
confined phase, namely the vison loop excitations of this 
phase are gapped and dilute. Although we do not yet 
have an explicit lattice model for this 3d SPT, the lattice 
model of this SPT phase must only contain terms that are 
even powers of bi and b[ in order to keep the U(l) symme- 
try: Hq = ^^ ■ —tb\ j6ij- + ■ ■ ■. Thus we can modify this 
Hamiltonian and couple 61 to a Z2 gauge field a^j defined 

on the links of the lattice: Hg = J2i j ~t^ij^i.i^i,j + ' ' • 
Now consider a long vison loop along z axis. This vi- 
son loop is bound with a half- vortex line of 61 (Fig. [ij, 
and the vison loop is the core of the half- vortex line. 
Along the vison loop (core of half- vortex line), since ri} 
and n? are zero, the effective Lagrangian along the vi- 
son loop only involves a three component unit vector 
n — (n^,n^,n^) ~ (Re[&2], Iin[62], </>)• The effective ac- 
tion along the vison loop reads 

Old = idltefn^dinf ^^, e,/ = l,2. (3) 

where I is the line coordinate along a large closed circle 
around the vison loop. 

In Eq. [31 Gid = TT is protected by time-reversal sym- 
metry. Under Zj, since a vortex of 61 transforms into 
an anti- vortex of 61, the derived Id 8— term changes its 
sign: 9id — ;■ — 0id, hence Eq. [3] is only time-reversal 
invariant at points Gid = irk with integer k. If we ig- 
nore the physical interpretation of the field n, this 1+ld 
NLSM at 9id = TT (Eq. [3]) can be used to describe the 
antiferromagnetic spin-1/2 chain, and based on the Lieb- 
Schultz-Mattis (LSM) theorem this l-|-ld system must be 
either gapless or degenerate [25|. When it is gapless, the 
vison loop is described by a 1 -I- Id conformal field theory; 
degenerate ground state can be induced by spontaneous 
time-reversal symmetry breaking along the vison loop. 

Notice that the vison loop is invariant under time- 



reversal transformation, because in one plaquette tt— flux 
and — TT— flux are equivalent. However, flux lines of other 
discrete gauge fields are not necessarily time-reversal in- 
variant, thus if we couple the same SPT phase to other 
discrete gauge fields, the line defects may be fully gapped 
without degeneracy. 

Example 2: 3d SPT with U{\) x Z2 symmetry 
This SPT phase is also described by Eq. [U with the 
following transformations of n: 



U{1) 



^ e'"b, 



,n 



n\ n" ^-n'',a^2,---5. (4) 



The C/(l) >^ Z2 symmetry is a subgroup of S0(5). Since 
6 — )► 6^ under Z2, the U(l) and Z2 symmetries do not 
commute with each other. 

The lattice model of this SPT phase can be constructed 
using bosonic rotor operator bj ^ exp(i0j) on lattice. 
The Z2 symmetry corresponds to the particle-hole trans- 



formation of b. 



and n° fields in the field theory 



correspond to the rotor density operator, which changes 
sign under Z2 particle-hole transformation. In order to 
keep the U(l) symmetry, the lattice Hamiltonian will 
only involve even powers of 6, thus we can couple b to 
a lattice Z2 gauge field. The rest of the analysis is very 
similar to the previous example: the half-vortex line of 
b bound with the vison loop will lead to a 1+ld 0(3) 
NLSM with Qid = tt along the vison loop, which must 
be either gapless or degenerate. Qid = tt is protected 
by the Z2 particle-hole symmetry: under Z2 transfor- 
mation Qid -^ —Old, because a vortex of b becomes an 
anti-vortex under particle-hole transformation. 
Exam^ple 3: 3d SPT with Z2 x Zj symmetry 
In this section we discuss line defects in 3d SPT phases 
with discrete symmetries only. Let us take Z2 x Z|" 
symmetry as an example. SPT phases with symmetry 
Z2 X Z2 have {^2)"^ classification according to Ref. [7|. 
These eight different phases can be built with three dif- 
ferent basic phases, one is the bosonic topological super- 
conductor with just Z2 symmetry. The other two corre- 
spond to the so called phase 1 and 2 of U{1) x Zj SPT 
phases in Ref. [13|, and by breaking the U(l) down to its 
subgroup Z2, the phase 1 and 2 in Ref. [l^ become SPT 
phases with Z2 x Z|" symmetry. All these phases are de- 
scribed by the same effective field theory as Eq. [1] with 
a different transformation of the 0(5) vector n under the 
symmetries. 

In this section we will take phase 1 of Z2 x Zj SPT 
phase as an example. In phase 1 of Z2 x Zj SPT phase, 
the vector n transforms as follows: 



Z, 



1 



n -> 



4,5. 



1. 



■5. 



(5) 



Presumably this SPT phase can be realized in a lattice 
spin system with a local Hamiltonian defined with spin 




FIG. 2: A vison line in a 3d SPT phase, (a), the two 2d 
boundaries exposed after cutting the SPT phase at 2: = 0. 
A{x) in Eq. [9] has a domain wall at 2: = 0. (6), the two 
boundaries can both be viewed as coupled Id wires that are 
described by 1+ld 0(4) NLSM with WZW term at level k = 
±1 (Eq.IIll). 



operators {S^ ^ S^ , S^) only. The Z2 symmetry can be 
viewed as the tt— rotation around 2;— axis. Based on the 
symmetry transformations, we can make connection be- 
tween field theory variables and lattice operators. For 
example, in phase 1 



n'^ix) ^ AaS^ 



Basy 



a = 4, 5; 



n''{x) ^ CbSJ + Db{S., xS,)-Sk + --- b ^ 1,2, 3(6) 

with real constant coefficients Aa, B\,, Cb and Db- 

Every term in the lattice Hamiltonian must only have 
even powers of S^ and S^ to protect the Z2 x Zj symme- 
try. Thus we can consistently couple S'^ and S"^ to a Z2 
gauge theory: Hg = E.,j -t<,^f ^} + ■ ■ ■ {a,b ^ x, y). 
With this coupling on the lattice, ■n!' with 6 = 4, 5 are 
coupled to the same Z2 gauge field 

The simple half- vortex line picture in the previous ex- 
amples is not totally applicable here, because there is no 
U(l) degree of freedom that can form a vortex around 
the vison loop. Thus let us instead consider the follow- 
ing structure: Cut the system open in the XY plane at 
z = 0, which will expose two boundaries (Fig. [2]3,). Both 
boundaries must have nontrivial spectrum, and they are 
both described by a 2-|-ld NLSM with 0(5) vector n plus 
a Wess-Zumino-Witten (WZW) term at level fc = ±1 re- 
spectively: 

Sa = d^xdr—ldafia) 
J 9 

± (fx du—eabcdenadxri'adyrf^dun'ldrrf^,{J) 



Here a = 1 , 2 denotes the top and bottom boundaries 
exposed. The 0(5) WZW term has level fc = 1 for top 
boundary (a = 1), and fc = — 1 for bottom bound- 
ary [a — 2) respectively. na{x,y,T,u) is an exten- 
sion of the space-time configuration na{x,y,T) that sat- 
isfies na{x,y,T,Q) ~ (0,0,0,0,1), and n„(a;,y,T, 1) = 
na{x,y,T). The boundary WZW term can be derived 
from the bulk O— term in Eq. [TJ because when 9 = 2n, 
the 3+ld bulk O-term can be written as 2+ld WZW 
terms at boundaries z — L and z = 0: O — term = 
WZWL,fe=i+WZWo,fc=-i. 

The symmetry of Eq. [7] needs to be reduced to the 
physical symmetry. Let us assume the system energeti- 
cally favors n^ over •n? , so we can integrate out nf and 112 
from Eq. [7]to obtain an effective action for 0(4) vectors 
^a = ("-ii^Qi^-Qi^-a)- If the system preserves the Z2 
symmetry, then the expectation values {n\) = (n|) = 0. 
Now after integrating out n^ , Eq. [7] is reduced to two 
0(4) NLSMs with a O-tcrm at 9 = ±tt: 



S = I d^x-{d^naf 
J 9 



± / d^x-—^eabcd£^LupKdf,n''^d^nldpn'^. (8) 

Here 9 = tt on the top boundary (or — tt on the bottom 
boundary) is protected by the Z2 symmetry. Detailed 
calculation of the 9— term at the boundary can be found 
inRef. [l3,|ll. 

Now let us reglue the two boundaries together, by turn- 
ing on the following coupling: 

3 
Sc = I d^xdr ^Bnl{x,T)nl{x,T) 



b=4,5 



A(x)n\(x,T)nl{x,T). 



(9) 



The coupling constant A{x) has a Id domain wall at 
X = z = 0: A{x) < for a; < 0, A{x) > for a; > 0. For 
the entire XY plane B <Q. This inter-boundary coupling 
corresponds to inserting a vison loop in the XY plane 
along the y— axis at a; = z = 0. For the half plane z — 
0, a; < 0, we can identify ni{x,T) = fi2{x,T) = n{x,T), 
and eventually the effective 2d action in this half-plane 
is an ordinary 0(4) NLSM with no 9— term. In the op- 
posite half-plane z = 0,x > 0, since A{x) > 0, we have 
{n\,n\,n\,n\) = {n\,n\,n\,—n\) = n, and the effective 
action for n in the half-plane a; > is an 0(4) NLSM 
with 9 = 27r: 



Sx>0 



I 



1 

:-| 
9 
i27r 



d^x-id^.n)'^ 



+ / d-^x—-^eabcdet,upn''d^,n''d^n''dpn'^. (10) 



Now the vison loop can be viewed as a Id domain wall 
of 9 between 9 = at a; < 0, and 9 = 27r at x > 0. 



Although both sides of the domain waU can be driven 
into a 2d gapped disordered phase without degeneracy, 
the domain waU must have nontrivial spectrum. Using 
the analysis in Ref . [2J| , if both sides of the domain wall 
are gapped, this domain wall (vison loop) is described by 
a 1+ld 0(4) NLSM with a WZW term at level-1: 



5 = / dydr -{d^,nf 



/[ i2Tr 

(fx / du-—^Tf-d^,-n!'dun''dpn'^eabcd<it,up{ll) 

It is well-known that this theory flows to a stable 1+ld 
SU(2)i conformal field theory fixed point under renor- 
malization group [26|, |27[, if the system has a full SO (4) 
symmetry. In our case, although the symmetry is much 
lower than SO (4), no linear term of n° is allowed by sym- 
metry. Any bilinear term of n° , even if it is relevant at 
the SU(2)i fixed point, will eventually lead to sponta- 
neous symmetry breaking and degenerate ground states. 
We seek for a more physical picture for the formal 
calculation above. Before regluing the boundaries to- 
gether, the boundaries are described by 0(4) NLSM with 
8 = ±7r (Eq. [S]). This theory can be viewed as coupled 
Id wires along y direction [28|, and every wire is a l-(-ld 
0(4) NLSM with a WZW term at level /fc = ±1 (Fig.[3>): 



-'a— 1,2, a;— J 



dydr -{d^fic,f± 



dydr / du eabcd^t.vpnl^d^j^Ti'^dyrf^dpn'l^li) 



m{-iy 

12772 _ 

If a direct inter-wire coupling X]q=i 2 S y "-(^ — J' J/j ''")« • 
n{x = j' + a, y, t)^ is turned on (a is the distance between 
nearest neighbor wires), each boundary reduces to the 
2-f'ld 0(4) NLSM with e = ±7r (Eq.© :2Sj. 

Now we glue the two boundaries together with a do- 
main wall of A{x). In the half plane cc < 0, since 
ni = n2 = n, two wires on top and bottom boundaries 
would trivially gap out due to their coupling with each 
other (their WZW terms cancel each other); however, 
on the other half plane a; > 0, due to the opposite sign 
of inter-boundary coupling, the WZW term of the top 
boundary wire x = j will cancel the WZW term of the 
bottom boundary wire x = j + a. Thus at the domain 
wall X — 0, there is one Id wire left which is not gapped 
by coupling with other wires. 

This picture is very analogous to coupling two spin- 1/2 
chains together. Let us consider two spin- 1/2 Heisenberg 
chains along x direction: H — J2a=i 2^1 ^j,aSj+i,a- 
At X < we couple the two chains antiferromagnetically: 
H' ~ J2j J'^j,i ■ "5^,2, while for a; > we couple the 

two chains ferromagnetically: H' = J2j ~J'^j,i ' Sj,2- 
Then the half-line a; < can be viewed as a trivial spin-0 
chain, while for a; > it is the Haldane phase of a spin-1 
chain. Then at the origin x = it is guaranteed to have 
a dangling spin- 1/2 doublet. 



The same kind of analysis and conclusion can be ap- 
plied to the phase 2 of SPT phase with Z2 x Zj sym- 
metry, where the 0(5) vector n transforms as Z2 : n^ ^ 
n\ n^ -^ -n!>,h = 2 • ■ • 5; Zj : n*^ -^ -n'^,a = 1---5. 
The only difference from phase 1 is that, in this case we 
need to couple rfi with 6 = 2 — 5 to the same Z2 gauge 
field. 

Example 4- Point defect in 2d SPT phase 
Let us now briefly discuss 2d SPT phases. A 2d SPT 
phase must have trivial spectrum in the bulk, but gap- 
less or degenerate spectrum on its Id boundary. But 
studies on quantum spin Hall insulator have suggested 
that if a point defect is created in a 2d SPT, this point 
defect might also change the spectrum. For example, if 
a quantum spin Hall insulator is coupled to a Z2 gauge 
field, then the vison excitation of this Z2 gauge field must 
carry a Kramers doublet 29|, |30| • 



Here we argue that similar effect also occurs gener- 
ally in 2d SPT phases. For instance, let us consider 2d 
bosonic SPT phase with C/(l) x: Z2 symmetry, which is a 
bosonic version of QSH insulator. This SPT phase is de- 
scribed by a 2-hld 0(4) NLSM with 9 = 27r T^ which 



involves a four component vector n — (n^. 



n ,n ,n 



')■ 



n^ +in^ is a boson rotor that transforms under U(l), and 
n'^,n^, n'^ all change sign under Zj. Let us couple n^ and 
n^ to a Z2 gauge field, and consider a vison at the origin 
of the 2d system. Then this vison is bound with a half- 
vortex of 6, which leads to a 0-1- Id 0(2) NLSM for n^ and 
n"* with Ood — TT a.t the origin: S = J dT^€ab'n°'drn'', 
a,b — 3, 4. This + Id model can be solved exactly, and 
its ground state is two fold degenerate, which is precisely 
a Kramers doublet. This degeneracy is again protected 
by time-reversal symmetry. Thus a vison excitation in a 
Z2— gauged bosonic quantum spin Hall insulator has the 
same behavior as the fermionic QSH state. 

If we break the U(l) symmetry down to Z2 (consider 2d 
SPT phase with Z2 x Zj symmetry) , we can still couple 
n^ and n^ to a Z2 gauge field. Now a lower dimensional 
version of Fig. [2] allows us to study the spectrum of the 
vison in the system, and the vison will still be two fold 
degenerate. The vison spectrum in this case can also 
be understood using the "decorated domain wall" con- 
struction of SPT phases discussed in Ref. [18| . In the 2d 
Z2 X Z2 SPT, a domain wall of the Z2 symmetry is a Id 
SPT phase with Zj symmetry, and after coupling the Z2 
part to a Z2 gauge field, a vison is the Od boundary of 
the Id SPT, thus it must be a 2-fold degenerate Kramers 
doublet. However, none of the defects in the previous 
cases discussed in this paper can be analyzed using the 
decorated domain wall construction. Our studies based 
on effective field theory are more general. 

In summary, we study the defects in SPT phases intro- 
duced by Z2 gauge field, and in all the cases discussed in 
this paper the defect (either line defect in 3d or point de- 
fect in 2d) has nontrivial spectrum. Our study not only 
reveals a new general property of SPT phases, it also 



suggests a possible way of classifying 3d Z^ topological 
order enriched by symmetry, based on the spectrum of 
its vison line. 
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